Resolving the inner active accretion disk 
around the Herbig Be star MWC 147 with 
VLTI/MIDI+ AMBER spectro-interferometry 
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Summary. We studied the geometry of the inner (AU-scale) circurristchar envi- 
ronment around the Herbig Be star MWC 147. Combining, for the first time, near- 

00 (NIR, K band) and mid-infrared (MIR, N band) spectro-interferometry on a Her- 

big star, our VLTI/MIDI and AMBER data constrain not only the geometry of the 
brightness distribution, but also the radial temperature distribution in the disk. For 

^ (-H our detailed modeling of the interferometric data and the spectral energy distri- 

^ I bution (SED), we employ 2-D radiation transfer simulations, showing that passive 

irradiated Keplerian dust disks can easily fit the SED, but predict much lower vis- 
ibilities than observed. Models of a Keplerian disk with emission from an optically 
thick inner gaseous accretion disk (inside the dust sublimation zone), however, yield 
^ a good fit of the SED and simultaneously reproduce the observed NIR and MIR 

1 I visibilities. We conclude that the NIR continuum emission from MWC 147 is domi- 

nated by accretion luminosity emerging from an optically thick inner gaseous disk, 
while the MIR emission also contains strong contributions from the outer dust disk. 

CO 1 Introduction and observations 

y—{ Herbig Ae/Be (HAeBe) stars are intermediate-mass objects which are still 

'O accreting material, probably via a circumstellar disk composed of gas and 

dust. Understanding the structure of these disks and the processes through 
which they interact with the central star is critical for our understanding of the 
^ formation process of stars. Since, until recently, the spatial scales of the inner 

k>( circumstellar environment (a few AU) were not accessible to infrared imaging 

observations, conclusions drawn on the 3-D geometry of the circumstellar 
material were mostly based on the modeling of the SED (e.g. Hillenbrand et 
al. [1]). However, these fits are known to be ambiguous (e.g. Men'shchikov 
& Henning [2]) and have to be complemented with spatial information, as 
provided by infrared interferomctry. 

The first systematic studies of HAeBes using the technique of infrared 
long-baseline interferometry have revealed that for most HAeBes the NIR 
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size correlates with the stellar luminosity L following a simple R oc L^/^ law. 
This suggests that the NIR continuum emission mainly traces hot dust at the 
inner sublimation radius. However, for more luminous Herbig Be stars, the 
NIR-emitting structure is more compact than predicted by the size-luminosity 
relation (Monnier et al. [3]). 

To further investigate the origin of the "undersized" Herbig Be star disks, 
we observed the Herbig Be star MWC 147 with the VLTI. First infrared in- 
terferometric observations of this star by Akeson et al. [4] indicated that the 
NIR-emitting region is surprisingly compact (^0.7 AU, assuming a uniform 
disk profile). In the course of three ESO open time programmes and using the 
8.2 m UT telescopes, we obtained seven VLTI/MIDI measurements and one 
VLTI/ AMBER measurement. The MIDI observations cover baseline lengths 
between 39 and 102 m. From the AMBER data, one wavelength-dependent 
visibility, corresponding to a 101 m baseline, could be extracted. For our mod- 
eling of MWC 147, we adopt the stellar parameters by Hernandez et al. [5], 
namely a spectral type of B6, a distance of 800 pc, a bolometric luminosity 
of 1550 Lq, a mass of 6.6 M©, and a stellar radius of 6.63 Rq. 

2 The power of joint NIR/MIR spectro-interferometry 

The VLTI instruments AMBER and MIDI combine the high spatial resolu- 
tion achievable with infrared interferometry with spectroscopic capabilities, 
measuring the fringe visibility as a function of wavelength. As circumstellar 
disks exhibit a temperature gradient, different spectral channels trace different 
spatial regions. Therefore, spectro-interferometric observations, which cover a 
sufficiently large wavelength range, can constrain not only the disk geometry, 
but also the radial temperature profile of the disk. For the investigation on 
YSO disks, the NIR and MIR wavelength regimes are particularly well suited 
since the NIR wavelength regime 2 /im) is most sensitive to the thermal 
emission from dust located at the dust sublimation radius {T w 1500 K, a 
few AU from the star), while MIR wavelengths 10 /xm) trace dust with 
a temperature of several hundred Kelvin, located a few 10 AU from the star 
(see Fig. 1). 




Fig. 1. Illustration of the inner environment of HAeBe stars. Due to the different 
temperatures, spectro-interferometry can disentangle multiple emission components. 



Resolving the Herbig Be star MWC 147 with VLTI spectro-interferometry 



3 



3 Radiative transfer modeling 

As a first step of analysis, we compared the interferometric data to commonly 
used analytic disk models with a simple temperature power-law {T{r) cx r~", 
with a = 3/4 or 1/2). Using the assumption that each disk annulus radiates 
as a blackbody, we can compute the wavelength-dependence of the disk size 
corresponding to these analytic models and find that they cannot reproduce 
the measured NIR and MIR-sizes simultaneously (see Fig. 2). 
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Fig. 2. Comparing the 
wavelength-dependent disk size of 
MWC 147 (as derived from one of 
our VLTI/AMBER+MIDI mea- 
surements) with the wavelength- 
dependent disk size predicted by 
standard temperature power-law 
disk models, we find that these 
models cannot reproduce our 
measurements. In this figure, the 
disk model was scaled to match 
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Wavelength [,mi the measured MIR size. 

Therefore, we applied a more sophisticated modeling approach using the 
mcsim_mpi 2-D radiative transfer code (Ohnaka et al. [6]). For each model, we 
first check the agreement with the SED of MWC 147 and then fit the spectro- 
interferometric visibilities (see Kraus et al. [7] for details). The dust density 
distribution of the accretion disk in our models resembles a flared, Keplerian- 
rotating disk with a radial density distribution of p{r) cx r^'^/^, which extends 
from the dust sublimation radius to 100 AU. In order to reproduce the shape 
of the SED, we found that, in addition to the disk, an extended envelope is 
required, for which we use the radial density distribution p{r) cx r~^/^. 

Fig. 3 shows model images, the SED, as well as the visibilities correspond- 
ing to our best-fit model of a passive, irradiated accretion disk. Although 
irradiated disks are able to reproduce the SED of MWC 147, they predict vis- 
ibilities which are much smaller than the measured visibilities {Xr ~ 26) and 
are therefore in strong conflict with our interferometric measurements. 

In passive circumstellar disks, the infrared emission is generally assumed 
to originate almost entirely from dust; the emissivity of the inner, dust-free 
gaseous part of the disk, at radii smaller than the dust sublimation radius, is 
negligible. In an actively accreting disk, on the other hand, viscous dissipation 
of energy in the inner dust-free gaseous part of the accretion disk can heat 
the gas to high temperatures and give rise to significant amounts of infrared 
emission from optically thick gas. The inner edge of this gas accretion disk is 
expected to be located a few stellar radii above the stellar surface, where the 
hot gas is thought to be channeled towards the star via magnetospheric accre- 
tion columns. While the magnetospheric accretion columns are too small to be 
resolved in our interferometric data (3 i?* correspond to 0.09 AU or 0.12 mas). 
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a) Image /f-band b) Image K-band c) Image 8^im d) Image 13^m 
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Fig. 3. Model images (a-d), SED (e), and NIR/MIR visibilities /-(^ corresponding to 
our best-fit radiative transfer image of an irradiated dust disk geometry. The poor 
agreement between measured and model visibilities (Xr ~ 26) indicates that passive 
disk models can be ruled out. 



infrared emission from hot gas between the dust subhmation radius and the 
stellar surface should be clearly distinguishable from the thermal emission of 
the dusty disk due to the different temperatures of these components and the 
resulting characteristic slope in the NIR- and MIR- visibilities. 

Since MWC 147 is a quite strong accretor (Afacc ~ 10~^ Mq yr"^; Hillen- 
brand et al. [1]), significant infrared emission from the inner gaseous accretion 
disk is expected. MuzeroUe et al. [8] found that even for smaller accretion rates, 
the gaseous inner accretion disk is several times thinner than the puffed-up in- 
ner dust disk wall and is optically thick (both in radial as well as in the vertical 
direction). In order to add the thermal emission from the inner gaseous disk 
to our radiative transfer models, we assume the radial temperature power-law 
by Pringle [9]. Including the accretion luminosity from an inner gaseous disk 
in the model strongly improves the agreement between model predictions and 
observed visibilities. With a flared disk geometry and an accretion rate of 
-^acc = 9 X 10~^ Af0yr^^, both the SED and the interferometric visibilities 
are reproduced reasonably well (x^ = 1.0, see Fig. 4). 



4 Conclusions and outlook 



Our VLTI interferometric observations of MWC 147 constrain, for the first 
time, the inner circumstellar environment around a Her big Be star over the 
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a) Image ff-band b) Image K-band c) Image 8^im d) Image 13^m 
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Fig. 4. Observables corresponding to our best-fit radiative transfer model with 
optically-thick inner gas disk (similar to Fig. 3), yielding good agreement (Xr ~ 1-0). 



wavelength range from 2 to 13 /im. We find evidence that the NIR emission 
of MWC 147 is dominated by the emission from optically-thick gas located 
inside the dust sublimation radius, while the MIR also contains contributions 
from the outer, irradiated dust disk. 

Our study demonstrates the power of infrared spectro-interferometry to 
probe the inner structure of the disks around young stars and to disentangle 
multiple emission components. Future investigations on YSO accretion disks 
will benefit substantially from the proposed 2nd generation VLTI instruments, 
such as MATISSE, increasing not only the number of recorded baselines, but 
also expanding the spectral coverage to the L and M bands. 
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